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Abstract 

A method is proposed for determining the response of the steady swirling potential flow 
past a blade row in an axial turbomachine to the most general type of i non-acoustic 5 
incident disturbance. The method is based on Goldstein’s (1978) decomposition of the 
disturbance velocity and only requires solving a linear inhomogeneous wave equation. 

It is believed that numerical solutions to this wave equation can be obtained more 
efficiently and will be more accurate than corresponding solutions to the linearized 
compressible Euler equations. 

1. Introduction 

The response of the steady compressible flow past a blade row in a turbomachine 
annulus to both steady and unsteady disturbances is of technological importance because of 
the effect such disturbances have on the performance of the turbomachine. The disturbances 
can arise from inlet-flow non-uniformities as well as blade-interaction effects and are often 
small enough to be described by linear theory. Furthermore, viscous effects are typically 
confined to thin boundary layers and wakes so the flow can be regarded as inviscid to good 
approximation. 

The behavior of small- amplitude disturbances in uniform flows is well understood 
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(Kovasznay 1953). The vorticity, entropy and pressure fluctuations associated with such 
disturbances are decoupled when the flow is unbounded and coupled through the blocking 
effect caused by the vanishing of the normal velocity on the impermeable surfaces when 
the flow is bounded. The situation is more complicated when the base flow is not uniform 
however. Goldstein (1978) presents a general theory for describing the behavior of small- 
amplitude vortical and entropic disturbances on arbitrary potential flows. This theory shows 
how the distortion of the vortical part of the disturbance velocity as it is convected by the 
steady base flow also gives rise to pressure fluctuations. In turbomachine flows, the vorticity 
and pressure disturbances are coupled through both the blocking and base-flow distortion 
effects. This is true even between the blade rows since the base-flow velocity there typically 
has a swirl component that is of the same order as the throughflow component. 

The purpose of this paper is to present a relatively simple method for determining 
the behavior of small-amplitude disturbances on the steady swirling flow past a blade row 
in an axial turbomachine. To this end, it is assumed that the base flow is irrotational, 
i.e. that the vortical and entropic parts of the motion are small enough to be completely 
accounted by the disturbance field. This assumption precludes the possibility of Rayleigh- 
type instabilities and allows one to use Goldstein’s (1978) decomposition of the disturbance 
velocity. This is done in §2 by writing the disturbance velocity as the sum of a vortical 
part u[ J ^ that is a known function of the imposed disturbance and a potential part 
that satisfies a linear inhomogeneous wave equation. Since the perturbation potential <j>\ 
determines the pressure fluctuations, it follows from the above discussion that 4>\ does not 
vanish in the upstream region where the affect of the blade row is negligible. In order to 
complete the formulation, an upstream boundary condition for <f > i must be derived. This 
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is done in §3 where an analytic solution for 4>\ is constructed for the case when there is no 
intervening blade row and the base flow is given by the superposition of a uniform stream 
and a line vortex. The results of §3 are then used together with the decomposition of the 
perturbation potential introduced by Atassi & Grzedzinski (1989) to obtain a problem for 
4>\ that is suitable for numerical integration. Finally, the results axe discussed in §5. 

2. Formulation 

Consider a steady inviscid non-heat-conducting compressible flow past a blade row in 
an axial turbomachine and suppose that a small- amplitude (steady or unsteady) disturbance 
is introduced in an upstream region where the effect of the blade row is negligible. It is 
assumed that the vortical and entropic parts of the motion in this upstream region are 
small enough to be completely accounted for by the disturbance field and further that any 
shocks waves that may exist always remain quite weak. The total flow will then be a linear 
perturbation about a steady irrotational base flow and wUl have velocity, pressure, density 


and entropy fields that expand as 


u = u o (r,0,z ) + u x (r,0,z,t) 4- 

(2.1) 

P = Po(r,M) + Pi(r,0,z,t) + ---, 

(2.2) 

P = po(r,0,z) + pi(r,9,z,t) + •••, 

(2.3) 

s = si(r,6, z,t) H , 

(2.4) 


respectively, where (r, 9 , z) axe cylindrical coordinates chosen so that the z-axis is coincident 
with the axis of the annulus and is pointing in the direction of the upstream throughflow, t 
is time and without loss of generality the constant base-flow entropy has been set equal to 
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zero. The origin of the (r, 0, z) coordinate system is taken to be in the plane of the blade 
row. The unit vectors in the radial, circumferential and axial directions are denoted e r , eg 
and e z , respectively, while the velocity components in these three directions are denoted 
u, v and w. Attention is restricted to a calorically perfect ideal gas so the specific heats 
c p and c v are constant and the temperature T is related to the pressure and density by 
P = (c P - c v )pT . 

The inner and outer radii of the turbomachine annulus can vary with z, but in the 
upstream region, it will be assumed that these radii take on the constant values r, and r 0 , 
respectively. Then the base flow becomes independent of z as z — * — 00 and is given by 

u Q = egToo/r + e z Woo, po/Po = constant, 
p^/cl = constant, Cq = (k - 1 ) ( H — uq-u q /2) , 

where 2;iToo and e z W 0 a are the constant circulation and throughflow velocity of the up- 
stream base flow, Co = (/cpo/po)* is the base-flow speed of sound, H is the constant stag- 
nation enthalpy of the base flow and n = c p /c v . 

The perturbation quantities in (2.1)-(2.4) are governed by the linearized compressible 
Euler equations which Goldstein (1978) shows can be reduced to a single linear inhomoge- 
neous wave equation when the base flow is irrotational. Using this reduction in the present 
formulation leads to 


Ui = u [ ^ + Vfc, 

(2.6) 

Do, 

Pl = ~ p0 Dt 

(2.7) 

_ J_Pi _ £1 
Pl ~ 4 Po Cp ’ 

(2.8) 

s i = — SaWoot), 

(2.9) 
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where D 0 /Dt = d/dt + tto'V is the convective derivative associated with the steady base 
flow, is the Kronecker delta tensor and the £; — t&Woot = (£ 1 , ^ 2 > ( 3 — Woo*) denote three 
functionally independent Lagrangian variables for the base flow. 

The vortical part of the velocity perturbation u[^ is given by 

u[ I] = Atfi - + ^ 5(6 - SisW^t), ( 2 . 10 ) 

where the repeated index j is summed from 1 to 3. The perturbation potential <f>\ satisfies 
the linear inhomogeneous wave equation, 

% (2 H) 

Dt \Cg Dtjp 0 po 

within the flow field and the inviscid boundary condition, 

n-V<h = — (2.12) 

on the impermeable boundaries of the flow where n is a unit normal vector. 

and £2 are two independent stream functions for the steady base flow and £ 3/^00 
is the drift function (Darwin 1954 5 Lighthill 1956) associated with that flow. Therefore the 
£i{r,8,z) are determined by 

uo -Vfc = SaWco, (2.13) 


subject to the upstream boundary conditions, 



(2.14) 


as z — * — 00 , where (2.14) is arrived at by using (2.5). These equations can be integrated 
to give the following analytic expression for £ 3 : 


6 = 2 + /_ 


Woo - w o (r a (£i,fr,0OA(fr>6»4 ,/ )> 2: »(£i>&> < ft / )) jjJ (0 I 5 j 

! 7"T T , f \ a / a > 1 + i/\N 19 > 5 \ * J 


00 K (r,(6, 6, <t>% Uii, 6, <f>% I 2 
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where 4>q is the base-flow velocity potential, i.e. uq = Vcfo, and r,, 9 a and z, give the 
variation of the cylindrical coordinates along a base-flow streamline and are determined by 
solving Ci = £i(r,0, 2 ), £2 = and <Ao = <fo(T,0,z) for r, 6 and z. 

The functions S and Aj of & - foWoo* are arbitrary apart from the restriction that 
they be periodic in & which follows from (2.14) and the 0-periodicity of the flow field. S is 
related to the entropy by (2.9) while the Aj can be related to the vorticity which is given 
as 

"> = vx '* 1 = { a ‘ + ■ (2 ' 16) 
Before deriving an explicit expression for Aj in terms of u>i, it is noted that the velocity 

decomposition (2.6) is not unique since an arbitrary function of £» - feWoo* can be added to 

<f>! and its gradient subtracted from to produce a new decomposition that also satisfies 

(2.7)-(2.12). In the present formulation, the velocity decomposition is made unique by 

requiring that the circumferential average of A\ and the circumferential derivative of A? 

both equal zero. Forming the scalar product of (2.16) and V£/, multiplying the result by 

the alternating unit tensor £/ mn , and using the identity £jhi£imn = 6jnf>km~f>jm$kn produces 


Wi At ~ Wk A ‘ “ n ‘ k ' 


(2.17) 


where the skew-symmetric tensor fij k is given by 


O (f _c w n ejkiuyVh 1 Wo, S') 

Qjk ^ 6i3W °°^ - V6 • V6 X V6 + 2c p d(tj,( k ) ’ 


(2.18) 


and is the Jacobian derivative. Then, in view of the restrictions on A\ and A 2 

given above, integration of (2.17) yields 

(2.19) 


Aj = J Qijd£\ + J (f?2 j — j) ^£ 2 ? 
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where a bar is used to denote circumferentially averaged quantities. Equations (2.9) and 
(2.19) show that the functions S and Aj can be completely determined from knowledge of 
the entropy and vorticity fields on a single surface of constant f 3 . 

The formulation is completed by specifying an upstream boundary condition for (2.11). 
As indicated above, the coupling of the vorticity and pressure fluctuations through the 
dipole-type source term /5 q 1 V -(poUj^ ) in (2.11) and the term in (2.12) continues 

in the upstream region. Consequently, <f> 1 does not vanish as z — ♦ —00. The appropriate 
upstream limit for 4 >\ is determined by first constructing a solution to (2.11) and (2.12) that 
applies when there is no blade row. 


3. Solution for <j> 1 in the absence of the blade row 


A formal analytic solution for the perturbation potential in an annulus with constant 
inner and outer radii where the base flow is given by (2.5) can be derived by considering a 
single harmonic component of the vorticity and entropy fluctuations, i.e. by letting 


(A 


( - \ 
Aitti) 


exp {i[fc 2 & + *3(6 - Woo*)]} 


Sit. 1 ) 


J 


(3.1) 


where £2 and fc 3 are constant wavenumbers (the former of which is limited to integer values 
due to the 0-periodicity of the flow field) and 


(3.2) 


A\ = 0 when = 0 

A 2 = 0 when ^ 0 

in order to satisfy the restrictions on A\ and Ai . Then since (2.11) and (2.12) are linear 
the solution for a general vortical/entropic disturbance can be obtained by superposition. 
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Since interest here is in the long-time response due to a small-amplitude disturbance 
forced at a fixed axial location say z = where — Zoo » 1, 4>i and d<f>\/dt axe taken to be 
zero at t = 0. Substituting (3.1) into (2.10) and the result along with (2.5) into (2.11) and 
(2.12) and taking the Laplace transform with respect to time then yields 

\ (iu, + t*o*V) 2 4>i - - V-(poV^) = + az \ (3.3) 

eg po W 

for r,- < t < t 0 and 

4-fa = —M^ + az ), (3.4) 

dr iv 

for r = r,, r a where 

&(r, e,z\u)= I™ e-' wt Mr, *, <)dt (3.5) 

JO 

is the temporal Laplace transform of <j>\, v = w + ^ 3^00 Is a shifted transform frequency, 


/ \ ^OO ^2 , L 

a(r) = -wZT>+ k * 


(3.6) 


is the r- dependent axial wavenumber of the vorticity and entropy fluctuations, and 

r oo 


it \ 1 f ® . yL Cso z 


rpo \dr 
+ i aA 3 + i k 3 


2c„ 


Woo r 3 

S 


) MO + ( 22 l - ihr ) (% Ai ) 


(3.7) 


and 

5 (r, z) = Ai + 2^^i 2 (3.8) 

are the harmonic coefficients of the dipole source term the radial com- 

ponent of the vortical velocity perturbation u[ ! \ respectively. 

In order to simplify the subsequent analysis, <j > i is decomposed as 



where the q n are determined by the inhomogeneous equations, 


\w 2 

/ ii/ . , d > 

) 2 -— — 

( 9 \ 

k 2 

d 2 ' 

oo 

r 2 

c 0 

VWoo ia + 

' rpo dr 

{ rp °¥rJ 

+ v>- 

' dz\ 


q n = tni /e 


1 az 


for Ti < r < r 0 and 


= ~^ aZ , 


(3.10) 


(3.11) 


for t = r,-,r 0 , and rp is determined by the homogeneous versions of (3.3) and (3.4). It 
follows from (2.7) that the q n account for the pressure fluctuations associated with the 
vortical entropic (or ‘non-acoustic’) part of the disturbance field while tp accounts for those 
fluctuations associated with the irrotational isentropic (or ‘acoustic’) part of that field. 

The solution for is found by introducing 

qi = f e ia ( r ') z qi(r,r',z\v)dr' (3.12) 

Jri 

into (3.10) and (3.11) to obtain 

+ 2i (a' - j; + (l - |^] 9i = -Kt ~ OM *)> ( 3 - 13 ) 


and 


respectively, where 


“ °- 


(3.14) 


T - 1 9 ( JL\-( — 

rp Q drV Po dr) \r 2 


+ o‘- 


2 Wl X 2 \ 

c 2 1 i 


, , » ^ m u 

A m,<r,u i( r ) — t 2 "t ° ' fy ’ 


(3.15) 

(3.16) 


a ' = a(r'), and 6(r) is the Dirac delta function. Then, in view of the linear dependence of 
/ on z, qi is written as 

<7i = Ki(r, r' \ v) + \K 2 (t, r' | /( r '> 2 )• (3-17) 
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By substituting this expression into (3.13) and (3.14) and using the results of appendix A, 
one can show that 


*i(r,r» = G(r,r'|^a>-W, (3.18) 

tf 2 (r,r» = ~G{r y r , \k 2 ,a',v-k 3 W 00 ), (3.19) 

where the Green’s function G(r,r'|m,a,w) associated with the operator L miCr , w is given by 
(A3). It is now a relatively simple matter to show, by using (A3) and (A 4), that 

q 2 = -e m ' z K\{r,r'\v) + iA^r, r'| i/) g{r',z) ^ , (3.20) 

since g is also linearly dependent on z. 

As indicated above, i/> arises from the ‘acoustic’ part of the disturbance field imposed 
at z = Zoq. The portion of the ‘acoustic’ field propagating inwards (in the direction of the 
throughflow) at z = Zoo can be specified arbitrarily, but the outward-propagating portion 
must be determined by the solution for No outward-propagating ‘acoustic’ waves arise in 
the absence of the blade row because the base flow is independent of z. Therefore y> - 0 for 
z < Zoo and the solution for can be found by taking the Fourier transform with respect 
to 6 and the Laplace transform with respect to z of the homogeneous versions of (3.3) and 


(3.4) to get 

L = -h for r, < r < r 0 

= 0 for r = Ti,T 0 

dr 

where 

^(r|m,<7,w)= f f e -1 t m ^ + CT (z _ z«O]^( r ,0,*|o;)dzd0, 

Jo J Zqo 

and 



(3.21) 


(3.22) 


h(r | m,<r,w) = 
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(3.23) 

Zoo 

is determined by the inward-propagating ‘acoustic’ field at z = z^,. The solution for ip is 
then just 

rj> = j G(r, r' \ m,(T,uj)h(r' \ m,<r,w)dr' (3.24) 

Jr, 

which when combined with (3.22) leads to 

■tp = — l — f ° ^ f e^ m ^" t ■ £^ ^ _^00 ^G(r,r / |m,<r,a;)fi(r'|m,<7,a;)do■dr , (3.25) 

(2 7r) J Ti m __ 00 Jc a 

where the integration contour C a must lie below the singularities of the integrand so that 
7 p = 0 for z < Zqo. The singularities correspond to singularities in G which, in view of (A 3), 
arise whenever 

D(m, <7,u>) = 0 (3.26) 

where D is given by (A 4). Therefore the singular points of G result from the discrete 
spectrum of the wave equation. 

It now follows from (3.5), (3.9), (3.17), and (3.20) that the solution for the perturbation 
potential is 

<*! = —/ e i[k20 + {v- Woo)*] J_£( r z | „)<!* + -L / e iw ^(r, 0, z | w)du; (3.27) 

2 tt J Cv iv 2nJc u 


W 2 
2i 00 


V>(r, 0, z | u;)d0 


where 

q = q\ + h = J e ia ' z K\{r, r' 1 1 /) + iK 2 (r, r ' 1 i/) /(/, z)dr' 

- e'“' 2 [ jr,(r, T>\v) + UTj(r, r< \ r) A] j(r', ,) , (3.28) 

■0 is given by (3.25), and in order for the solution to be causal, i.e. 4>\ = 0 for t < 0, the 
integration contours C v and C u must lie below the singularities of q/ u and V>, respectively. 
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The long-time behavior of the ‘acoustic’ part of the perturbation potential is determined 
by those singularities associated with the discrete spectrum that lie on or below the real- 
v axis. The discrete spectrum associated with the base-flow solution (2.5) is considered 
in Kerrebrock (1977) and the interested reader is referred to that paper for a detailed 
discussion. However it should be noted that the ‘unstable’ modes found in Kerrebrock 
(1977) are excluded from (3.27) by causality. For the remainder of the present study, 
attention will be restricted to ‘non-acoustic’ disturbances so will be set equal to zero. In 
addition to contributions from the discrete-spectrum singularities, the long-time behavior 
of the ‘non-acoustic’ part of the perturbation potential includes a contribution from the 
singularity at v — 0 which determines the response at the forcing frequency of the imposed 
vortical/entropic disturbance. 

4. Solution for 4n in presence of the blade row 

The solution for the perturbation potential derived in the preceding section provides 
the upstream boundary condition needed to complete the formulation begun in §2. Once 
the base-flow solution (tto. Po> Po) 1S known and the & have been determined from (2.13)— 
(2.15), the response of an actual swirling turbomachine flow to the most general type of 
‘non-acoustic’ incident disturbance can be found from (2.6)— (2.12). However difficulties 
arise when one tries to solve (2.11) and (2.12) numerically subject to (3.27) as z — * — oo. 

It is clear from (2.15) that £3 becomes singular at the base-flow stagnation points, 
i.e. the points where uq = 0. The zeros of the Jacobian determinant <9(£i , £ 2 , 4>o)ld(r,6, z) 
also give rise to singularities in £3 but these points can be shown to coincide with the zeros of 
uq. The singular behavior of & is, in general, transferred through u[ T) to the perturbation 
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potential fa thereby making (2.11), (2.12) and (3.27) unsuitable for numerical integration. 
Atassi & Grzedzinski (1989) overcome this difficulty by expressing the solution for fa as 

fa = <f\ (6 - £.3 Woof) + 4>X(r, e , z,f), (4.1) 

where <pi is a known function that is singular at the base-flow stagnation points and <f>* is 
a non-singular function that satisfies a linear inhomogeneous wave equation. 

For simplicity, it is assumed that the base flow has only one stagnation point which is 
denoted P . The relevant results for flows with multiple stagnation points and/or stagnation 
lines can be obtained by superposition. The surface formed by the base-flow streamlines 
issuing from P is denoted E and, in general, consists of a portion E\> that lies on an 
impermeable boundary of the flow and a portion E w that lies within the flow itself. By 
assuming E is smooth near P, Atassi & Grzedzinski (1989) show that will contain all 
the singular behavior in <f>\ if 

n-(t4 J) + V<p!) -> 0 and -u 0 -(«{ 7) + Vy?i) — (4.2) 

ic v 

as E is approached and furthermore, that it is sufficient to take 

¥>! = - f 3 ~ Wcot {^(fi 0) ,4 0) ,0 + exp [o(& -r- Woof)] 1,6,0} d£\ (4.3) 

where a = W^{d\uo\/ dn)p, 

B(ti - ^,-3 Woof) = - a [(6 - *x 0, )Ai(fc - 6* Woof) + (6 - ~ *3 Woof)] 

+ ~ ^Woo f ) - A 3 (rf 0) , 4 °\ 6 - Woof), (4-4) 

and a superscript (0) indicates evaluation on S. 
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It now follows from (2.6), (2.7) and (4.1) that the velocity and pressure perturbations 
produced by the imposed disturbance can be expressed as 

ui = Ji* + V<£*, (4.5) 

and 

Pl = -Pojfitl, ( 4 - 6 ) 

where D 0 /Dt = d/dt + u 0 -V is the convective derivative associated with the steady base 
flow. The known solution for tq } is 

u (R) = Aj(ti - + |r 5 (6 - feW’oo*) + » ( 4 - 7 ) 

where (2.5), (2.9), (2.14) and (2.19) can be used to relate the functions S and Aj to the 
vorticity and entropy fluctuations of the incident disturbance at z = z^. 

Substituting (4.1) into (2.11) and (2.12) and using (4.2) one finds that 

yr - T-V.(poV^) = y Vip 0 u[ R) ), (4.8) 

Dt\clDt J po Po 

which must be solved subject to 

n>V<£* = 0 on Sb, 

A(n-V^) = 0 on S w , y ( 4 - 9 ) 

ra-V^f = -n-u^ on 77, 

where A denotes the jump across E w and 77 denotes the parts of the impermeable boundaries 
not included in E b . The solution for 4>* must also satisfy the upstream condition 

<f>t -» ~<fii + 4><x> (4.10) 

as z _► -oo, where (f>oo is given by the long-time behavior of the right-hand side of (3.27) 
but with $ set equal to zero. 
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5* Discussion 


The determination of the response of a steady swirling potential flow past a turboma- 
chine blade row to the most general type of ‘non-acoustic’ incident disturbance has been 
reduced to solving the linear inhomogeneous wave equation (4.8) subject to (4.9) and (4.10). 
This reduction required modifying the upstream boundary conditions given in Goldstein 
(1978) because the vorticity induced pressure fluctuations do not vanish in the upstream 
region of the present problem. These fluctuations result from the coupling of the vorticity 
and pressure fields through the blocking and base-flow distortion effects and persist in the 
upstre am region due to the vanishing of the radial velocity on r — , r 0 as well as the 

swirling nature of the base flow. 

A numerical solution to (4.8) — (4.10) for the flow past a blade row is beyond the scope 
of the present study but some of the results obtained from the solution that applies in 
the absence of the blade row will be discussed. The results of §3 imply that the long- 
time behavior of the ‘acoustic’ part of the disturbance field is given by a sum over the 
discrete spectrum associated with the wave equation but that the long-time behavior of the 
‘non-acoustic’ disturbance field can not be given by such a sum. 

Tan & Greitzer (1986) considered the behavior of steady non-axisymmetric distur- 
bances on compressible swirling flow in turbomachine annuli. The authors found that the 
am plitudes of these disturbances decrease as the axial distance from the point of excitation 
increases. The large-z behavior of q is given by (B 6) of appendix B for the case when k 2 ^ 0 
and this result together with (3.27) shows that the amplitudes of the non-axisymmetric ‘non- 
acoustic’ disturbances decrease like z~ l a s z —* oo. In contrast, (3.7), (3.8) and (3.28) imply 
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that the amplitudes of the axisymmetric ‘non- acoustic’ disturbances grow linearly with z 
in this limit. 

The difference between the large- z behaviors of axisymmetric and non-axisymmetric 
parts of the ‘non-acoustic’ disturbance disturbance is a consequence of the way in which 
the relative motion between particles frozen in the base flow (i.e. purely convected with 
it) affects the solutions to (2.11) and (2.12). Equation (2.14) shows that the base-flow 
particle trajectories in the upstream region are helical and that the particle velocity has 
an axial component equal to Woo and a circumferential component inversely proportional 
to the (constant) radius of the trajectory. Therefore the relative motion only occurs in the 
circumferential direction and then only between particles of differing radial positions. This 
means that the terms in (2.11) and (2.12) that involve derivatives with respect to r become 
dominant as z becomes large but only when the ‘non-acoustic’ part of <f > i varies with 0. 

Appendix A. The Green’s function G 

In this appendix, a Green’s function G associated with the operator L is con- 
structed. The governing equations for G(r, r' | m, a ,u) are taken to be 

L m ,<T,uG = —6(r - r') (A 1) 

for Ti <r <r 0 and 

jr-G = 0 (A 2) 

dr 

for r = r;,r 0 . Let G\(r \ m,<r,u) and G 2 (r \ be two linearly independent homoge- 

neous solutions to (A 1) that satisfy the initial conditions Gi = 1/rpo, dG\/dr — 0 at r = r, - 
and G 2 = 1/rpo, dG 2 jdr = 0 at r = r 0 . Then the solution to (A 1) and (A 2) can be written 
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as 


G(r, r' | m,cr,u)= - * ^ ^ t G\ (r< | m, <r, u ))G 2 {t > \ m, <r, u) (A 3) 

v ' ' D(m,cr,u) 

where 

D = --^Gi(r 0 |m,<7,u;) = -Q^G 2 (ri 1 m,a,u;) (A 4) 

is the product of rpo and the Wronskian of G\ and G 2 and r< (r>) is the lesser (greater) 
of r and r\ 


Appendix B. The large- z behavior of q 


In this appendix, the behavior of q in the limit of z — > oo is determined for the case 
when £2 7^ 0. Considering first, (3.12) is rewritten as 


4i = jT e ia ' z gidr' + jP e^'^dr' 


(B 1) 


so that qi(r,r',z\i/) can be continuously differentiated with respect to r'. Integrating by 
parts and using (3.2), (3.7), and (3.17) then shows that 

J e ia ' z q\dr' = - J~/ a ' Z { /(r', *)iSTi(r, r» - M{r')K 2 (r,r' \ v)(3' 

- [ii(r , )A'i(r,r , |i/)]| +0(z~ 2 ) (B2) 


as oo where 


3( r \ - — - o r °° k - 
P{r) ~ dr “ Voor®’ 


(B3) 


and ft = P(r'). Equations (3.18), (3.19) and (A3) imply 


K x {r,r'\v) 


r'=r+ 


= tf 2 (r,r» 


r'=r+ 


= 0 , -^iKi{r,T'\v) 


= 1 , 


r'=r+ 
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and 


^7#i(r,r» 


J’ - { v^)t [ ’' Vtl(r ' )1 * ,(r ’ r ' 1 ^ - p ' K2(t ’ r<l L_* ’ 


so substituting (B2) into (B 1) yields 


x * _iaz ; , 1 „ia'z 

<7l = TS“ e Al + ’^aT e 

i/9 z i/9'z 


i/3'zAiir') + ia'A 3 (r') 




S'i(r.r'k) 


r 0 


+ 0(0 (B 4) 


asz-4W. Then, since (3.20) reduces to 

h = -e IQ,2 Ai(r , )ii'i(r,r , |i/) 


(B5) 


when ^ / 0, 

* = -A-J az ^ + ^-/ a ' z 


i /9z 


. w . 

ia'A 3 (r') + iA: 3 — —S(r') 


*i(r,r» 


+ 0(z~ 2 ) (B 6) 


as z — » oo. 
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